Preisach images of a simple mechanical system 
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This work is an an early stage of a larger project aiming at answering the question whether or 
not the Preisach map is really fingerprinting magnetic materials. More precisely, we are interested 
whether Preisach model of magnetic hysteresis indeed contains any physics or is just a convenient 
modeling tool. To this extent we study a very simple mechanical system, thus fully tractable, sub- 
jected to the external force. Despite of its simplicity, our model captures all the fundamental features 
of real magnetic materials, namely their hysteretic behavior, coercivity, remanent magnetization and 
saturation at high fields. Both the overall shape of major hysteresis loop as well as First Order Re- 
versal Curves (FORC's) are reproduced quite correctly; they are very similar to those observed in 
magnetic materials. The model essentially consists of a single, spring loaded, rigid and rotative bar 
with non-zero friction torque. The length of a projection of this bar onto the direction of an external 
force is identified with magnetization. The friction torque and the spring constant are the only freely 
adjustable parameters of our model. Here we investigate, and present, their influence on the inferred 
Preisach maps. 
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I. MOTIVATION 

It took more than ten years to put into experimental 
practice jl| the idea of the so called FORC's (First Order 
Reversal Curves) first introduced by Mayergoyz Q in 
1986. Since then the FORC diagrams became more or 
less standard method of analysis of magnetic materials 
exhibiting hysteretic behavior. The FORC diagrams are 
nearly identical with Preisach Q maps, except for the 
boundary part (necessarily equal to exactly zero in FORC 
picture) , corresponding to the so called reversible part of 
magnetization curve. For this reason we will be using 
both terms as if they were equivalent. 

In the beginning, the shape of Preisach map was es- 
sentially guessed. According to some theoretical hints 
given by Neel (1949-1952), such a map should consist 
of a single peak. The parameters of this peak, assumed 
most often as being Gaussian, Lorentzian or lognormal, 
were later refined in such a way as to reproduce the ex- 
perimental data. The FORC technique, in contrast, does 
not assume any particular shape of a Preisach map and 
therefore reproduces empirical facts more accurately. It 
is really good news for engineers, but physicists rather 
would be interested what is the meaning of various details 
observed in result of FORC-type analysis. More specifi- 
cally, it is very important to be able to draw some defi- 
nite conclusions concerning various magnetic interactions 
dominating in a given material. 

The published data with recovered FORC diagrams 
are quite ubiquitous but lack systematic attitude. The 
experimental data are sometimes accidental (rock sam- 
ples) , or only a very limited set of samples is available. In 
addition, the experiment itself is rather time-consuming. 



On the other hand, simulated data are easy to generate 
and, what is even more important, the analyst 'knows ev- 
erything' about her artificial 'sample' in advance. Manip- 
ulating various parameters is straightforward and prac- 
tically unrestricted. 

In this situation, it seems natural to develop systematic 
studies in that matter, starting from the simplest mod- 
els and making them more realistic (and more complex) 
later. Here we present our first attempt in this direction. 



II. THE MODEL 

In our model (see Fig. [1} a stiff rod can be rotated 
around the axis located at (0,0). The rod is connected 
with the rotation axis by means of a spring. An exter- 
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FIG. 1. A sketch of the mechanism imitating magnetic sys- 
tem. In the absence of friction, the equilibrium position of 
a rotary rod (thick line) is at 9 — 0. The external force F is 
always parallel to the a;-axis. 

nal force is applied either parallel or antiparallel to the 
x-axis, so the rod's deflection 9 from vertical direction 
may never exceed w/2. A static friction between the rod 
and its rotation axis is present. No movement is possible 
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whenever the torque produced by spring and the external 
force has lower magnitude than the static friction torque, 
i.e. when 

|r s + r F K|r f |, (1) 

where the spring torque r s = — k9, with the dimensionless 
parameter k being called spring constant from now on, 
and the external force produced torque is rp = (p.xF) z = 
fi sin(6> — |) = (i cos6>, where /«(= 1) is the length of 
a rod, and, finally, T{ is the static friction torque. 

When the relation ([T]) is not satisfied, the rod rotates 
in the direction dictated by the sign of r s + tf until the 
equilibrium condition 

\t s + t f \ - \ T {\ =0 (2) 

is reached. 

In our mechanical model of magnetic hysteresis the ex- 
ternal force F is identified with the exciting field H , and 
the projection of a rod onto the direction of external field, 
equal to /i sin 9, corresponds to (reduced) magnetization, 
since fi=l. 

ILL THE RESULTS 

We present three selected hysteresis curves in Fig. [5] 
The shape of FORCs, perfectly horizontal, may be a lit- 
tle surprising at first sight, since our FORCs differ signifi- 
cantly from the majority of published data. Yet, contrary 
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FIG. 2. Hysteresis in the system under study. Parame- 
ters (spring constant, friction) are: (70,70) for dotted loop, 
(150,80) for dashed one, and (180,50) for the loop with 99 
first order reversal curves shown. 

to intuition, our simulations are not completely unreal- 
istic. One might expect horizontal reversal curves only 
for very hard magnetic materials, but the curves simi- 
lar to ours have been recently observed for arrays of Ni 
microwires Q electrodeposited in nanoporous alumina 
templates. The overall shape of our recovered FORC di- 
agrams is roughly the same for sensible combinations 
of adjustable parameters. Three examples are shown in 
Fig. El 




FIG. 3. FORC diagrams for the same combinations of spring 
constant and friction as in Fig. [2] From top: (70,70), 
(150,80), and (180,50). Note: the z-coordinate is scaled to 
keep the height of peaks fixed, but noisy contour z — is 
omitted to show more clearly the peak's shape. 



It is worth noticing that the peak is not a convex body 
- note the presence of 'arms' and 'folded ellipses' as con- 
tour lines. This fact alone invalidates some quasi one- 
dimensional, simplified models of the Preisach maps. It 
is also in sharp contrast with other analytical model of 
magnetic hysteresis, proposed by Takacs Q and shown 
in Fig. |31 The results of more quantitative analysis are 
shown in following figures. 

The position of the highest point on a map, more pre- 
cisely its coercivity field coordinate, is strictly propor- 
tional to the friction but is insensitive to the spring stiff- 
ness (Fig. [SJ . The error bars there are taken as the field 
values nearest to the crossing of major hysteresis loop 
and M = axis. 
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FIG. 4. Preisach image for Takacs model of magnetic hys- 
teresis p|. Note the very regular peak shape. 



The peak's amplitude, on the other hand, is a decreas- 
ing function of both parameters. It seems inversely pro- 
portional to the spring constant (see Fig. [6j) . The depen- 
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FIG. 5. The coercivity is strictly proportional to friction. 

dence on friction is less clear, due to inaccuracies gen- 
erated by imperfect smoothing procedure, resulting in 
a slightly warped peak shape, as one can see in Fig. [3] 
Nevertheless, it is quite likely that the peak's height is 
inversely proportional to the friction as well. 



IV. CONCLUSIONS 

Our results show that various approximations of the 
peak shape, existing in literature, are probably too sim- 
plistic or even incorrect. 

We occasionally observe small negative regions on ob- 
tained Preisach maps, especially at their edges, near co- 
ercivity coordinate equal to zero (see the top of Fig. [3]). 
We are strongly convinced that it is an artifact, due 
to the numerical differentiation and/or smoothing pro- 
cedure, rather than the real thing. 
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FIG. 6. Main peak's amplitude vs. spring constant for dif- 
ferent friction constants. Solid line shows l/(spring constant) 
dependence. Friction is (from top): 10, 20, 30, 40, 50, 60, 70, 
80, 100, 120, 140, and 160. The lines connecting points are 
eye- guides only. 

Negative regions are often observed both in real and 
simulated @ magnetic systems. This feature is most 
likely related to the multimodal character of a potential 
energy landscape in real systems (pinning sites, defects, 
domain structure), while in our model the potential en- 
ergy is unimodal. 
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